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Abstract 

This paper deals with the stabihty properties of a closed market, where capital and labour 
force are acting like a predator-prey system in population-dynamics. The spatial movement 
of the capital and labour force are taken into account by cross-diffusion effect. First, we are 
showing two possible ways for modeling this system in only one country's market (applying a 
simple functional response and a HoUing-type ratio-dependent response as well), examining 
the conditions of their stability properties. We extend the ratio-dependent model into two 
countries common market where two kind of cross-diffusion effects are present, and find those 
additional conditions, whose are necessary for the stability of the global common market 
besides the stability of each countries local markets. Our four-dimensional model highlights 
that a hectic movement of the capital toward labour force can cause a Turing instability. 



Key words and phrases: Capital-labour system; functional response; ratio-dependence; cross- 
diffusion; Turing instability 
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1 Model set up and preliminaries 



It is known by Goodwin [T3] that some economic situation can be described by a Lotka-Volterra 
predator-prey system. Capital plays the role of prey species while labour force is the predator 
species. Instead of a simple Lotka-Volterra system it is often used a logistic growth rate of capital 
|13] . The inhomogeneous distribution of capital and labour force has an important role as it 
was shown in [11], [1]. The model (1.3) of [11] was modified by cross-diffusion and this effect 
was modeled by a PDE. In [3] the same model was studies in a patchy environment, namely, 
the diffusion was modeled by an ODE. It was shown in both cases that cross-diffusion usually 
causes Turing instability. In our paper we state a different model in which both kind of cross- 
diffusion effect are present. Basically we are interested in such cases when the diffusion cannot 
cause instability in our higher dimensional model. 

In papers [11], [1] the kinetic system was modeled by a predator-prey situation as well. We 
revise the kinetic system according to the following situation. Let us consider a country, where 
a firm takes place. It has incentives to expand using its capital stock to invest. For these new 
firms labour force is required. The total labour force, i.e. the number of those employed and 
those unemployed, is denoted by v{t) at time t. The company supplies free jobs, which number 
is denoted by u{t) at time t and it is proportional to the invested capital, thus, both labour force 
and free jobs are attracts each other for sustenance. We assume that the quantity of free jobs 
would grow according to the logistic law if there was no labour force available. The labour force 
decreases the per capita growth rate of free jobs proportional to the quantity of labour force. The 
factor of proportionality is given by m > 0. This is the rate by which the labour force is filling in 
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the free jobs. This leads to the following differential equation: 



u{t) 



u{t)r 1 



u{t) 



K 



) 



mu{t)v{t) 



(1) 



where r is the natural per capita growth rate of free jobs and K > is its theoretical eventual 
maximum. Equation ([1]) is exactly the same as (1.1) in [TT]. We assume that the per capita growth 
rate of labour force is increasing proportionally with respect to the free jobs and we consider the 
case when the factor of proportionality is equal to m. Moreover we suppose that in the absence 
of free jobs, the labour without revenue is unable to maintain itself, leading to its declension. The 
death rate of labour force is denoted by > 0. This leads to the differential equation: 



This equation is different from the equation (1.2) given in [TT]. In population dynamics the function 
mv{t) is called functional response. It is often fails to be linear. Later we also show such model in 
which the functional response is more complicated. But in this simple case the following system 
of differential equations characterizes the situation: 



Comparing model ([3]) to the one given in [TT], we can see that both model is a predator-prey 
situation where capital is represented by prey species while labour force is by predator species. 
The only difference between them that in the model given in [TT], both prey and predator can 
survive without the other, while in ([3]) the predator dies out in absence of prey. Moreover, model 
([3]) is the same as model (2.5) without delay given in [13]. As in [TT] we modify model ([3]) by taking 
into account that the economy of the country is not concentrated in a point but is disturbed in a 
spatial bounded domain. It is assumed that capital and labour force is moving around freely in 



v(t) = mu{t)v{t) — dv{t). 



(2) 




(3) 
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the country and that the capital investment and the conditions of hving are the same everywhere. 
Besides, we assume that the economy is closed, meaning that there is no in- and outflow of capital 
and labour force at the boundary. In order to be able to perform explicit calculations we consider 
the spatial domain, abstractly, as one dimensional and identify it with the interval [0, L] C M, 
L > 0. From now on, u{t,x) and v{t,x) represent the densities of free jobs and labour forces, 
respectively, at time t > at the point x G [0, L] of the domain. We assume that the movement of 
capital and labour force follow Pick's law, i.e. both substances flow away from places where their 
densities are high, towards places where their densities are low. Furthermore, both substances are 
influenced by each other: larger quantity of free jobs attracts labour force rather then smaller, and 
new firms are more likely to be founded in places, where the density of labour force is higher. The 
velocities and the direction of flowing can be represented by diffusion coefficients and their signs, 
contained by the diffusion matrix: 



where du and d22 denote the self-diffusion of capital and labour force respectively, du the velocity 
of capital-flowing towards labour force, and d2i shows the largeness of the migration of labour 
force towards free jobs and dij > 0, i, j = 1,2, denote the self-diffusion effect (e.g. [12] )• In order 
to obtain a well posed problem, necessary to have det D > 0. (The same diffusion effect was 
considered in [11].) Now we can set up our first reaction-diffusion model: 




(4) 



Utit,x) 



u(t, x)r ( 1 



u{t,x) 
K 



) 



mu{t, x)v{t, x) + duu"^{t, x) — di2v"^{t, x) 



(5) 



v't{t,x) 



mu{t, x)v(t, x) — dv{t, x) — d2iu"x(t, x) + d22v'^^(t, x) 



<(^,0) = <(t,L) = = <(t,0) 



<(t,L) 



(6) 
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where boundary condition ([6]) presents the closeness of the market. The initial conditions are given 



It will be shown that at certain values of dij, i, j = 1,2 a Turing bifurcation may occur as it often 
happens in similar situation. 

This paper is organized as follows. In section |2] we are giving an overview of the stability 
behaviour of model ([5]). It will be shown that this model has some disadvantageous properties. In 
section [3] model ([5]) will be modified by a so-called ratio-dependent functional response and the 
effect of the changing of the diffusion coefficients will be studied to the stability behaviour. In 
section m we are extending the basic model to a patchy environment illustrating the outer movement 
of capital and labour force between two countries without and with inner migration. We state and 
study a four-dimensional PDE model which is capable to describe the inner and outer movements 
of the market at the same time. The model ( IT3l) in section [3] is a well known model and model 
( 122|) in section |4] is also known. We repeat here partly the known stability investigation of them 
and we emphasize some stability conditions by other aspects in order to understand better our 
new model (150]) in section SJ In section \5\ we are summarizing the main conclusions of the study. 

2 Stability of the simple market 

The kinetic system of ([5]) has three equilibria, namely: 



by 



u{0,x) = U{x), v{0,x) = V{x), xe[0,L], U,V E £^{[0, L],R). 



(7) 



So = (0,0); Ei = {K,Oy, E={u,v) 



where 




dr + Kmr 



Km? 



8 
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and V > if K > —. Each of them can be examined, however we have special interest in E 
(given by ([8])), because this is the only one which is not placed on the boundary of the plane. The 
equilibrium of the kinetic system is also a constant solution of (|5]). The interaction matrix of the 
kinetic system linearized at E is: 



Thus, the following theorem holds: 

Theorem 2.1. Suppose, that condition ^ is satisfied; then the equilibrium E = {u,v) of the 
kinetic system of ^ is asymptotically stable. 

Theorem 12.11 shows that system ([5]) is a too simple model because if a positive equilibrium 
exists then it is necessarily stable. (We note here that if condition ([9]) holds then matrix Jg is 
signstable e.g. p!5].) 

Now we are linearizing the reaction-diffusion system ([5]) around the constant solution {u,v), 
using the notation: 




It is easy to see, that this is in fact stable, if 



K > —. 

m 



(9) 



pit, x) = u(t, x) — u; 



q{t, x) = v{t, x) — V. 



We get 



p't{t,x) 
q't{t,x) 



{r-^ - p{t, x) - {mu)q{t, x) + dup'^^{t, x) - di2q"^{t, x) 
{mv)p{t, x) + {mu — d)q{t, x) — d2ip'^x{t, x) + rf229L(^i ^) 




(10) 



The boundary conditions are obtained: 



V.it^ 0) = pl(t, L) = = gl(t, 0) = gl(t, L). 



(11) 



Cross-diffusion modeling in macro-economical model 



8 



Solving ( ITOl) by Fourier-method, it is easy to see (e. g. ^0\, [H]) that the stabihty behaviour 
depends on the stabihty of the foUowing matrices: 



A - A.D 



— Xkdu 



Km 



-d + Xkdi2 

Matrix A — A^D is stable if the following conditions are satisfied: 

dr 



, A, = ( — A: = 0,1, 2, 



Tr(A - DAfc) = Tr(A) - X^TrCD) 



Krn^ 

>o 



Xk{dn + d22) < 



>o 



det (A - DAfc) = det A +A? det D -A^f--— ^22 - ddii + r 1 - -— 1 ^12) > 0. 
— V — ■' — V — ' Km \ Km 

>o >o ^ 



Since det A > because of ([9]) and we assumed initially, that det D > 0, the sufficient condition 
to get a positive determinant for all Xk is 

dr I d \ 

d22 - dd2\ + r 1 - -— dx2 < 0. 



Km 



Km 



(This is the case when det (A — DA^) is a stable polynomial.) Expressing di2 from the above 
condition, we obtain: 



(ii2 < 



dr 
Km 



d22 + dd2i 



\^ Km) 



(12) 



Thus, we have arrived at: 



Theorem 2.2. // conditions ^ and fT^) hold then the constant solution (0, 0) of the linear 
problem [Id\) -fTl\ ) is asymptotically stable. 

Remark 2.1. According to [8] (or to [18]), the asymptotic stability of the constant solution [u, v) 
of system ([5])-([6]) is implied. 
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Increasing the diffusion coefficient du, exceeding the critical value: 



d 




■21 



12, 



\ Km J 



where condition f ll2p does not hold any more, equilibrium (m, v) of system ([S])^® may undergo a 
Turing-bifurcation [19] and became unstable. More precisely, if di2 > '^i2crit then there exists 
for which det (A — DA^) < 0. Examining the sign of the diffusion coefficients in condition (1121) . it 
is clear that higher values of ^21 and ^22 are supporting the stability of the system, while di2 has 
a destabilizing effect as it increases. Hence, our calculations lead us to the following conclusion: 
rapid movement of labour force towards free jobs and fast moving of labour force to areas where 
its densities are low from places where its densities are high are advantageous for the stability 
of the market, while the hectic movement of capital stock (free jobs) towards places where the 
labour force density is higher acts negatively, even has the ability to make the market unstable. 
Interesting fact that diffusion coefficient du does not have an influence in stability in this case. 
Namely, the movement of free jobs from the places where its density is high towards the places 
where its density is low has no importance. This latter is also a disadvantageous characteristic of 
this model. Thus, we turn to study our model using another functional response. 

3 Model with Holling— type ratio— dependent response 

In population dynamics a Holling-type ratio-dependent functional response is used in that cases 
when the competition for food is very sharp, see for example in [1], [6], [7], [9], [5], [IE], [H]. It 
is reasonable to apply this kind of functional response when the labour force need to search for 
free jobs. According to this, the following model is relevant when the competition for free jobs is 
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intensive: 

u't{t, x) = u{t, x)r (l - ^) - + duKAt^ ^) - di2v"xit^ ^) 

v[{t,x) = m^^0i^-dv{t,x)-d2iu'^^{t,x)+d22v'^^{t,x) 

with the boundary conditions (Q and the initial conditions ([7]). Parameter a > is called half- 
saturation constant (e. g. [16]) and m > is the maximum growth rate of labour force. Of course 
the model is relevant if the maximal growth rate of the labour force is larger than its death rate, 
namely 

m-d>0. (14) 




The equilibria of the kinetic system of (fT3|) are: 

Ei = {K,0); Er = {u,v), 

where 

, .Kid — m + ar) Kim — d){d — m + ar) . 

[u, v) = ( , — ) (15) 

ar a^ar 

Remark 3.1. Equilibrium Eq = (0,0) is missing, because certain denominators would be equal 
to zero. This problem could be solved by continuous extension of the model (e. g. [16]), however 
we are interested in E,.. 

Necessary to have the equilibrium E^ positive, which requires the following condition: 

r > . (16) 

a 

This means that the natural per capita growth rate of free jobs (capital) has to be higher than the 
ratio of the maximal growth rate and half saturation constant of the labour force. This is natural 
because high growth rate of labour force (relatively to the half saturation constant) would result 
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a too fast falling in free jobs. 

Following the method given in section [21 linearizing ( 1T3|1 in Er (given by f fTSj) ) we obtain the 
interaction matrix in the following form: 



ma m 
(d—m)'^ d{m—d) 



This is stable if det Jg^ > and Tr{J^J < 0, namely 

d(m — d)(d — m + ar) 

det Jg = > 0, 

am 

which holds if and only if conditions f lT6|) and f fT^ hold. 

^ , (a—l)d'^ — adm + m(m — ar) (a — l)d'^ — adm + md 

TriJsJ = ^ < 

am ^v^rl am 

by ^6} 

<0 by IT^ 

(a — — c?m) (a — l)d {d — m) 



< 



am am 

holds if 

a > 1. (17) 

In population dynamics the low value of parameter a refers to a predator species which tries to 
ensure its survival by having a relatively low growth rate and consuming less. It is a so-called 
K-strategist (e.g. [10]). In our case the meaning of ffTTj) will be that we apart from low willingness 
to work, low desire to work, low interest to work or discharged worker. Summing up the conditions 
we have obtained: 

Theorem 3.1. Suppose that conditions [14\ ), and (11) hold; thus equilibrium Er (given by 

fT5\) ) of the kinetic system of ( fi^) is asymptotically stable. 
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Now, just as in section Ej we solve the linearized model of (IT3l)-(l6l) by Fourier-method. We 
get that the stability behaviour of the spatially constant solution depends on the stability of the 
following matrices: 



A, - DAt 



-d2 



r - Xud 



fe"ii 



m n 1^ 



fc"22 



A, = A^ fc = 0,l,2, 



In order to be this matrix stable the following conditions need to be satisfied for all Xk- 

Tr{Ar - DAfc) < 0; det (Ar - DA^) > 0. 

From condition (fT7|) we know that Tr(Ar) < 0, using this: 

rr(Ar - DAfc) = Tr(Ar) - A,, [du + ^22). < 0. 

>o 

Let us consider 



det (Ar — DAfc) = det Ar + A^ det D — Afc(am(i22 + ar22'^ii + '^ri2'^2i + ar2i'^i2) > 0, 

where Ar = [ory], detD > 0, and if conditions ffT^ and ffTB]) hold, we have det Ar > as well, 
thus, a sufficient condition to have a positive determinant is 

d22 + (2^22 '^11 ~l~ '^ri2'^21 + Q,r2i di2 < 0. (18) 

In this case det (Ar — DA^) is a stable polynomial. Substituting the given entries of Ar and 
expressing ^12: 

m + d 



, ad ad^ 
"12 < ;«ii + 7 ;t^"2i + 



+ 



amr 



m — d {m — dY \ m — d (m — dY 

Applying again the theorem of Casten-HoUand [8] we get: 



d^ 



22- 



(19) 
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Theorem 3.2. Suppose that conditions [l4\ ), ( Tf^) - ( i7), ( fJ^) hold then the constant solution Er 
of system / l7^) -/f^) zs asymptotically stable. 

Remark 3.2. The coefficient of d22 takes positive values in fll9p if 



amr 
{m — d) 



>m + d, (20) 



or 



m — d ,^ d . , 

r> 1 + — 21 

a m 

holds. From this the positivity of follows. In this case < 0. This means that we are 
outside of the so-called Allee-effect zone of free jobs. (See the meaning of this zone in [TU], [TB] in 
connection with a prey species.) Outside of the Allee-effect zone the number of free jobs are not 
too low. 



Increasing diffusion coefficient di2, exceeding a critical value, where condition ( 1191) does not 
hold any more: 

^ '^d ^ ad^ ^ / d + m amr 



system f lT3l) -fl6|) may undergo a Turing-bifurcation and further increase of di2 causes instability. 
(Similarly as in section |2l) Examining the diffusion coefficients in condition f fT9|) . it is obvious 
that higher values of du, ^22 and d2i (if ( l20l) holds), are supporting the stability of the market, 
while an increase in coefficient di2 affects negatively on it. Namely, fast movement of ffims towards 
rarely industried areas supports the stability of the market, so as rapid migration of labour force 
towards places where its densities are low in case of greater unemployment and where free jobs are 
abundant. Nevertheless, hectic flow of the capital towards places where labour force is dense acts 
negatively and also able to violate the stability of the system. 
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4 The patchy model 

In this section we study the cross-diffusion effect between two countries. They will be modeled 
by patches e.g. [21 [31 H]. We apply these methods for our cases as follows: Let us consider two 
countries as two patches, providing the following: the inner migration remains still inside each 
patch, but a new movement of capital and labour force has taken into account between the two 
countries. We are curious about the conditions of the stability of the common market, providing 
that each market is also stable, and about how the diffusion affects on the obtained stability. The 
Holling-type ratio-dependent model introduced in section [3] will be used. 

4.1 Without inner migration 

First we consider the patchy model as two points, thus we neglect the inner migration of the two 
substances. The model can be written as follows: 

u{l,t) = f{u{l,t),v{l,t)) + 6i{pi{v{2,t))u{2,t)-pM^,t))u{l,t)) 

v{l,t) = g{u{l,t),v{l,t)) + 62{p2{u{2,t))v{2,t)- p2{u{l,t))v{l,t)) 

u{2,t) = f{u{2,t),v{2,t)) + 6,{pM^,t))u{l,t)-pM^,t))u{2,t)) 

v{2,t) = g{u{2,t),v{2,t)) + 62{p2{u{l,t))v{l,t)-p2{u{2,t))v{2,t)) 

where numbers 1, 2 in the arguments refer to the number of patch, thus, those mean the number of 
the country, while / and g are the same functions as in the kinetic system of the ratio-dependent 
model f lT3|) . namely: 

f(u,v) = ru(l-^)-m^ 
q(u,v) = m-^ dv 

Pii P2 £ are the migration functions. Similar model was studied in |1] but now the reaction term 
(l23l) is different. These latter functions describe the diffusion of capital and labour respectively 



(22) 



(23) 
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between the two countries (patches), 6i is the diffusion coefficient, denotes the velocity of this 
movement. Functions pi,p2 > 0, and p'i,P2 < 0- (Here p[{v{i,t)) = ^^^^^^wT^' P^i.^i.h'^)) = 



dp2{u{i,t)) 
du{i,t) ' 



dv(i,t) 

1,2.) These positive decreasing functions expressing the facts that the migration of 



the labour is (free jobs are) higher into a country from the other if the number of the free jobs are 
(labour is) less in the other country. These functions model the cross-diffusion effect between the 
countries. We say that the cross-diffusion is strong if \p^\ is large. If pj = 1 then we have mere 
self-diffusion. Since problem (|22l) has the same kinetic system as problem ( fT3l) -(j6l) if we restrict 
it to a single country, all the equilibria and stability conditions will remain unchanged. Let us 
denote the interaction matrix of system (|22|) at Ep = {u, v, u, v) by Ap (where u, v are given by 
ffTSD). we get: 



A. 



v 



drll Or-12 

flr21 0^22 

ttrll arl2 

a21 ar22 



( 



J 



— r 



v 



am, 






m 
(m.-d) 












-d^ 



{d-mf 






_^ 

m 

d{m—d) 



\ 



J 



Theorem 4.1. Suppose that conditions [I4\ ), fTB) . hold (for both countries), thus, the equi- 
librium Ej. is asymptotically stable in each markets, which implies the asymptotic stability of equi- 
librium Ep = {u,v,u,v) of the kinetic system of ^2E) (where {u,v) is given by l[T5\) ). 

Now denote the coefficient-matrix of the linearized reaction-diffusion system (l22l) at Ep by 



(Ap + r), where 



^ -6ipi{v) -Sip[{v)u Sipi{v) Sip[{v)u ^ 

-^2P2(^)^ -^2P2(m) hp'2{u)v S2P2{u) 

Sipi{v) Sip[{v)u -Sipiiv) -dip[{v)u 

y S2p2{u)v 52p2{u) -52p'2{u)v -52p2{u)- J 
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The characteristic polynomial can be determined by row- and column operations: 
det (Ap + r - yul) 



5ipi{v) 



r22 



62p2{u)v 



- 62P2{u) 

6ip[{v)u 

S2P2{u) 



cl+c3; c2+c4 



^22 



6ipi{v) 6ip[{v)u 

P S2P2{U)V hp2{u) 

Orii - hPl{v) - P - Sip[{v)u 

Oral - S2P2iu)v 0^33 - S2P2{u) - P 

6ipi{v) Sip[{v)u 



S2p'2{u)v 



S2P2{u) 



"rii 



a. 



ri2 



«rii - Sipi{v) - p a,.i2 - 6ip[{v)u 



a. 



r2l 



'•T22 



r3— rl; ri—r2 



a. 



Til A* '-^ri2 






P ar2i — S2p'2iu)v ar22 — 52P2{u) — p 

5ipi{v) 5ip[{v)u 

52p'2{u)v 52P2{u) 

r-ii - 25ipi{v) - p a^ia - 25ip'^{v)u 



det 



a^ai - 252p'2{u)v 022 " '252P2{u) - p 

o-ni ~ fJ- 0^12 o.rii — 2(5ipi(v) — p ar^2 '~ '26ip[{v)u 

det 

ar2i a.r22 — P Cir2i — 2(52^2 ('0)^ ^22 " 252^2 (^i) " P 

The first determinant in the previous product is the characteristic polynomial of the kinetic system 
of f[T^ . thus, according to theorem 13.11 it is stable if conditions dH]), ( 1TB]) and f lT7|) are satisfied. 
Consequently we have to determine the cases where the second determinant, meaning a polynomial, 
will be stable. Let us denote the matrix hy B = whose characteristic polynomial is given by 



det 



— 25ipi{v) — p a^ij ~ 2Sip[{v)u 
ar2i - 252P2iu)v ar22 - 252P2{u) - p 



(24) 
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If conditions 1^, (EH), p- < 0, > 0, i = 1, 2 and 

ttr,, - 25ip[{v)u < (25) 

are satisfied tlien B is sign-stable e.g. [U] . Substituting and expressing 5i from we 
obtain: 

5i < -f- . (26) 

2p[[Vjum 

Parameter 6i shows the velocity of the capital (free jobs) movement between the countries. Ac- 
cording to fl26l) . 6i has to be under a limit, namely too fast flow of capital is not allowed. With 
this consideration, we have arrived at: 

Theorem 4.2. Suppose that conditions (I^, (T^, (n^), ^2^, p- < 0, Pi > 0, z = 1,2 hold for 
each patch and ( f^) also holds for system ( fll|) then equilibrium Ep = {u,v,u,v) of system 
(given by fT5\)) is asymptotically stable. 

This theorem highlights, that for the stability of the common market, relatively slow foreign 
firm-movement and a significant number of free jobs are required besides the original stability 
conditions for each market individually. Consequently, the hectic cross-country flow of the capital 
has to be avoided again. 

If this particular capital investment is fast enough i.e.: 

2p\[y)um 

then the common market may lose its stability. 

The stability of fl24l) can be ensured by the following way, as well. The entries in the main 
diagonal in B are negative, thus, the determinant of B has to be positive in order to get stability: 
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det B = detAr - 2{5ipi{v)ar22 + 52P2iu)arii - 5ip[{v)uar2^ - 52p'2iu)vari^) + 4:5i52{pi{v) p2{u) - 
p'i{v)up2{u)v) . If (JH]) and f l2T|) are satisfied, Ar is signstable, tlius, det^l^ > 0. If 



1 > pi^^pmi (27) 



then {pi{v)p2{u) — p[{v)up2{u)v) > 0. If 



< ^ (28) 



then )a.r22 + 52P2(^)arii — )Mar2i "~^2P2('^)^^ri2) < 0- After short calculations we obtain 

( 128|) in the following form: 

- i < (29) 
Now we can summarize our result in the following theorem: 

Theorem 4.3. Suppose that conditions ^2^, H^i, ^29^, p\ < 0, pi > 0, i = 1,2 hold 

for each patch then equilibrium Ep = {u,v,u,v) of system ^2E) (given by ^^^) is asymptotically 
stable. 

Example 4.1. In the case of pi{v) = fg, p2{u) = ^ a,/3 > 1 conditions ([27]), ([29]) hold 
automatically. 



4.2 With inner migration 

Now imagine model (122|) with the inner migration of each patch, provided that the economy is 
closed. The spatial domain inside a country is one-dimensional again, as it was in case of the 
model f|T3|) . The "length" of the z-th country is denoted by Li, and [0, Lj] C M, Lj > 0, z = 1,2. 
We assume that in each country, capital and labour are increased and decreased by the in- and 
outflowed quantity. Besides there is inner migration as well. (The situation is similar to the 
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one which occurs often in hydrodynamics when a source and a "sink" are in a closed domain.) 
Moreover, we suppose that in case of migration an employee prefer such area in the other country 
which is situated the same distance from the border of the countries. More precisely, let us consider 
an employee living in place x G [0, Li]. If it will move to another country then it will choose exactly 
the place y = Xj^ E [0, L2] first, to find job. Considering these ideas, the complete model can be 
written as: 

u'^{l,t,x) = dn{lK,{l,t,x) - d^2{lK,{l,t,x) + f,{u{l,t,x)^ 

v'^{l, t, x) = -d2i(l)<;.(l, t, x) + d22(l)<.(l, oc) + g,{u{l, t, x), v{l,t, x), u{2, t, x^^),v{2, t, f^x)) 

u[{2, t, y) = rfii(2)<;(2, t, y) - du{2)v'^A2, t, y) + /2(n(l, t, T/f^), v{l, t, y^f^),u{2, t, y), 0(2, t, y)) 

v[i2, t, y) = -rf2i(2)«^;.(2, t, y) + ^22(2)01(2, t, y) + ^2(^(1, t, y|^), t;(l, t, 2/g), ^(2, t, y), 0(2, t, y)) 

(30) 

where 

/i(ix(l, t, x), v{l,t, x), n(2, t, xf^), 0(2, t, xf^)) 

= /(u(l, t, x), ^;(1, t, x)) + 5i(pi(t) (2, f, xf^))n(2, t, xf^) - pi(^;(l, x))u{l, t, x)) 
^i(w(l, t, x),vil, t, x), {i(2, t, ), 0(2, t, )) 

= 5'(w(l,i,a;),z;(l,t,x)) + ^2(^2(0(2, t, xf^))?)(2, t, xf^) - ^2(^(1, a;))v(l, t, x)) ^^^^ 
/2(ix(l, t, y^^),v{l, t, 2/f^), 0(2, t, ?/), 0(2, t, y)) 

= /(0(2, t, y), 0(2, 2/)) + 6i{pi{v{l, t, yj^^)Hl,t, y^) - pi(0(2, t, y))0(2, t, y)) 
hHl.t, y^J, v{l,t, ?/f^), 0(2, t, y), 0(2, t, ?/)) 

= ^(0(2, t, y), 0(2, t, y)) + 52(P2(m(1, t, yj^^))v{l,t, y^J - p2{u{2, t, y))v{2, t, y)), 
where f{u{l,t,x),v{l,t,x))) and f{u{2,t,y),v{2,t,y))) are given by ( !23l) . The boundary condi- 
tions are similar to (E]): 

<(1, t, 0) = <(1, t, Li) = = <(1, t, 0) = <(1, t, LO, ^^^^ 
0' (2, t, 0) = 0' (2, t, L2) = = 01(2, t, 0) = 01(2, t, L^). 
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Let us change variable y by Xj^ and introduce u{2,t,y) = u(2,t,Xj^) = u{2,t,x), v{2,t,y) 
v{2,t,xj^) = v{2,t,x), dij{2) = dij{2){j^y. Thus, we obtain: 

u[{l,t, x) = dii{l)u'!^^{l, t, x) — di2{l)v'^^{l, t, x) + t, x), t, x), ti(2, t, x), t;(2, t, x)) 

z;^(l,t,a;) = -(i2i(l)tii'a;(l, t, x) + d22{^)vl^{l,t,x) + ^1(^(1, t, x), t, x), m(2, t, x), 'u(2, t, x)) 

?/;(2, t, x) = dn(2)<,(2, t, x) - (ii2(2)<,(2, t, x) + JiHl.t, x),v{l,t, x),u{2, t, x),v{2, t, x)) 

v',{2, t, x) = -d2i(2)<,(2, t, x) + d22(2)<',(2, t, x) + ^2(^^(1, t, x),v{l, t, x), m(2, t, x),v{2, t, x)) 

(33) 

where 

fi{u{l,t,x),v{l,t,x),u{2,t,x),v{2,t,x)) 

= f{u{l,t,x),v{l,t,x)) + 5i(pi(f(2,i,x))n(2,t,x) - pi{v{l,t, x))u{l,t, x)) 
gi{u{l,t,x),v{l,t,x),u{2,t,x),v{2,t,x)) 

= g{u{l, t, x), v{l,t, x)) + 52{p2{u{2, t, x))v{2, t, x) - p2{u{l, t, x))v{l, t,x)) 

(34) 

/2(m(1, t, x), t, x), ti(2, t, x),v{2, t, x)) 

= /(m(2, t, x), ^;(2, t, x)) + Si{pi{v{l, t, x))u{l, t, x) - pi{v{2, t, x))m(2, t, x)) 
^2(^(1, x), v{l,t, x), u(2, t, x), i'(2, t, x)) 

= 5'(m(2, t, x), t;(2, t, x)) + 62{p2{u{l, t, x))v{l, t, x) - ^2(^(2, t, x))v{2, t, x)), 
where f{u{i,t,x),v{i,t,x))) are given by (!23|) (2 = 1,2). The boundary conditions are similar to 

m- 

t, 0) = t, Li) = = t, 0) = Li), 

(35) 

<.(2, t, 0) = <(2, t, Li) = = t4(2, t, 0) = <.(2, t, Li). 
This latter is equivalant to u'y{2, t, 0) = Uy{2, f, L2) = = ^^^(2, t, 0) = £"^(2, L2). The non-trivial 
equilibrium point Ep of system (l22l) is a spatially constant solution of system (l33l) (and of (l30l) . as 
well) , because 
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(where {u,v) is given by ( fT5l) ) and substituting Ep = {u,v,u,v) into ( IMll : 
fi{u,v,u,v) = f{u,v)+5i{pi{v)u- pi{v)u) =0 



=0 =0 

gi{u,v,u,v) = g{u,v) +52 {p2{u)v - p2{u)v)^ = 

=0 



>, ^ = 1,2. 



=0 



Let we introduce the following notations: U 



u{l,t,x), v{l,t,x), u{2,t,x), v{2,t,x) 



and F{U) 



fl, 9l, f2, 92 



Thus, we can transform model ([3SD to the following form: 



Ut = D,U,, + F{U) 
I 



(36) 



where the diffusion matrix is: D4 



V 



rfii(l) -rfi2(l) 

-rf2l(l) ^22(1) 

dii(2) -rfi2(2) 

-41(2) 6/22(2) j 



Linearizing sys- 



tem (136|) at Ep and introducing the notation V = U 



we obtain: 



y; = D4\/: + (Ap + r)\/. 



(37) 



We can solve equation fl371) by Fourier-method again. It is easy to see that the stability of Ey 
depends on the stability properties of the matrix (Ap+F — AfeD4), where = (f^)^' A; = 0, 1, 2, ... . 
In order to get its characteristic polinomial in a much simplier form, we use the following row- and 
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column operations: 



det (Ap + r - AfcD - fil) : 

- <5ipi(ii) - Afedii(l) - /.t ari2 - hp'i{v)u + Xkdi2{l) 5ipi{v) 5ip[{v)u 

ar2i - S2p2{u)v + Xi.d2l{l) ^22 - <52P2(«) - Afcd22(l) - M &2p'2{u)^ &2P2{u) 

Sipiiy) &ip'-^^{v)u aril - 5ipi(ii) - Afe(iii(2) - ^ - 5ip'i {i))M + Afcdi2 (2) 

S2p2(u)v &2p2{u) arji - 52P2 + Afed21 (2) ar22 - 52P2{«) - Afed22{2) - ft 

aril - Afec(ii(l) - ari2 + -'*fciii2(l) (SipiCii) 5ipi{j;)ji 

«r2i + Afed2l(l) ar22 - -'*fc'^22{l) - A* feP2(")^' <52P2(«) 

-Afc(dii(2) -cin(l)) -Afc{di2(2)-di2(l)) - 25ipi{j;) - Afedii(2) - ai2 - 25ip; + Afedia (2) 

cl + c3; -Afc{d2i(2) -ci2i(l)) -Afe{d22(2) - ^22(1)) 0,-21 - 252P2{«)« + Afed2i (2) ^^22 - 252P2{«) - Afcd22{2) - ft 
c2 + c4 

r3 — rl; 

r4 — t2 



Let us consider the special case when there is no difference between the velocities of the diffusion- 
coeffitients inside the two countries (patches). This is valid for the self- and the cross-diffusion. 
Thus, 

d,^{2) = dij{l) := d,j, 1,3 = 1,2. (38) 
In this case the latter determinant is equal to: 



det 



Urii — Afedll — p ari2 + AfcCil2 

Oral + Afed2i ai-22 " -^fe<i22 — P 



det 



aril - 25ipi{-i;) - X^dn - p - 25ip[{v)u + Xkdi2 

ar2i - 252P2{u)v + Xkd2i ar22 - 252P2{u) - Afed22 - P 



(39) 



Taking into account the signs of the coefficients: dij > 0, pi > 0, 6i > 0, p'liy) <Q, p^iu) and 
assuming f l2T|) we obtain: 



det 



Ajirfii — p ari2 + Aj.(ii2 



<0 

Oral + Afed21 ar22 — Afed22 — P 

>0 <0 



det 



Orii - 2Sipi{v) - Afedii -p ai2 - 25ip{{v)u + X^di 



arai - 252p'2{u)v + Xkd21 Qraa — 252p2{u) — Xkd22 — P 
>0 <0 



In order to get a stable polinomial, these two determinants have to be separately stable polinomials 
depending on p. Since the entries in the main diagonal are negative, the two determinants (without 
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subtracting fx in the main diagonals) have to be positive. Namely the following two conditions 
have to be satisfied: 



(40) 



(aril ~ ^kdll){ar22 — Afc(i22) — («ri2 + ^fcf^l2) (cir2i + ^kd2l) 

= Xl det D - Afc(arii(i22 + CLr^idn + 0^2^21 + ar2idi2) + det Ar > 0. 

(a^ii - 2dipi{v) - Afcdii)(a,.22 - 262p2{u) - Xkd22) 
-(ai2 - 26ip[{v)u + Afc(ii2)(ar2i - 262P2{u)v + \kd21) 

= Xl det D 

— Afc [(0^11^22 + (^r22 

dll + 0^12(^21 + «r2i'^12)] (41) 

+2[52P2{u)du + 6ipi{v)d22 + S2P2{u)vdi2 + Sip'i{v)ud2i]) 
+ {detAr - 2{5ipi{v)ar22 + S2P2{u)ar^^ - 6ip[{v)uar2^ - 52p'2{u)var^^) 

+4.5i52{pi{y)p2{u) - p\{v)up'^{u)v)) > 0. 

In connection with (140|) it can be seen that if (fT9|) is satisfied then it is a stable polynomial, thus 
it is positive for all A^ > 0. 

In connection with (gl]) it can be seen that if (119]), ([27]), (EH]) and 

52P2{u)dii + 6ipi{v)d22 + 52P2{u)vdi2 + 5ip[{v)ud2i > (42) 

are satisfied then it is also a stable polynomial, thus it is positive for all A^ > again. 
Condition f H2]) can be rewritten in the following form: 

S2P2{u)du + 6ipi{v)d22 + Sip[{v)ud2i , . 

"12 < r ■ l^dj 

-62p2{u)V 

We have arrived at: 
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Theorem 4.4. Suppose that conditions (T^, {SP, (E^j, g^, p- < 



0, Pi > 0, 2 = 1,2 hold for each patch then equilibrium Ep = {u,v,u,v) of system ^3^) (given by 
( [73]) ) is asymptotically stable. 

This theorem means that if we consider a country as a patch it is required not to have a too low 
number of free jobs in a country (according to fl2Tl) ): the hectic flow of the capital towards places 
where labour force is dense has to be avoided inside the country (according to ( |T9i) and ( l43l) ): the 
migration functions pi, i = 1,2 has to be characterized by ( (!27l) . (l29l) ) that describe the inflow of 
capital and labour into the two countries (patches), and those give a limit for inflow of capital, the 
two countries have similar market practices (according to 0381) ) then the common market is stable. 

Comparing the conditions of Theorem 14.31 and Theorem 14.41 it can be seen that it is not enough 
to guarantee the stability of the markets between the two countries. If we consider the spatial 
movements inside the countries then stability conditions are required more stronger. In other 
words, not properly designed or supervised trade mechanism between the respective countries 
may as well cause the destabilization of each market, irrespectively of their own stability. Now 
the reversed case might also happen, namely due to a change in any (or both) market designs, 
the originally balanced (stable) joint market can lose its stability. All this means, that internal 
(external) market changes cannot occur without the distortion of the external (internal) market. 

Remark 4.1. By Casten-Holland (1977) [8], the asymptotic stability of the constant solution 
(m, V, u, v) of system (130|) is implied. 

Increasing diffusion coefficient di2 system ( !30l) may lose its stability. 
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5 Conclusion 

In section |3] the stability of the market in a country was considered in the case of diffusion of 
capital and labour. It was established that the hectic movement of the capital stock toward places 
where labour force is dense acts negatively and able to destabilize the system. We determined 
a critical value for the cross diffusion coefficient (^12) under which the market is stable in the 
country. Then in section H] two similar country was considered with their own stable markets. We 
studied the effect of the movement of capital and labour between the countries to the stability 
behaviour. In this case the hectic movement of the capital inside a country has to avoid even 
more. It may be necessary to reduce the inner flow of capital toward labour more significantly. 
This result is not surprising, since capital flow, namely a new investment, in order to obtain new 
jobs has to be slower than the movement of the labour. As a conclusion, the common market can 
be unstable if the inner movement of capital is sufficiently fast. This movement has the same effect 
on the individual market of each countrie, hence the destabilization of each countrie may affect the 
destabilization of the common market, and vice versa: instability of the common market makes 
the individual markets to be unstable. The wise moderation of capital is useful in all cases. Our 
results are applicable in case of two similar countries, in those where the velocity of capital flow 
and labour movement is approximately equal in the two countries, namely where there is similar 
market behaviour. Thus, our investigation can be applied for EU members, but not for an EU 
member and China common market. In this latter case we have to reject condition ( l38l) because 
the market practices are different and the stability investigation is become more complicated. 
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